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We examine the appearance of a spontaneous bulk spin current in a triplet superconductor in 
contact with a metallic ferromagnet. The spin current results from the spin-flip of Cooper pairs 
upon reflection from the interface with the ferromagnet, and is shown to display strong similarities to 
the spontaneous charge current in a Josephson junction. We express the spin current in terms of the 
Andreev reflection coefficients, which are derived by the construction of the quasiclassical scattering 
wavefunctions. The dependence of the spin current upon a number of parameters is investigated, in 
particular the orientation of the magnetic moment of the ferromagnet, the exchange splitting, the 
temperature, and the orbital pairing state of the triplet superconductor. 
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I. INTRODUCTION 

The interface between a singlet superconductor (SSC) 
and a ferromagnet (FM ) is an ideal setting to explore the 
antagonistic relationship between these two phases*^ It 
is now more than 25 years since the first theoretical in- 
vestigations,— but the physics of SSC-FM interfaces con- 
tinues to fascinate and surprise J^A^L^l^lL^JMl A 
key feature of such systems is the existence of an un- 
conventional proximity effect ^^ i 10 ' 11 ' 13 ' 14 in contrast to 
the spin-singlet pairing state of the bulk SSC, the spin- 
splitting of the Fermi surface induces spin-triplet correla- 
tions in the FM. Although this effect occurs at all metal- 
lic FM interfaces, the character of the induced triplet 
pairing correlations is determined by the strength of the 
exchange-splitting in the FM ^ 10 ' 13 whether the system 
is in the ballistic or the diffusive limits ^ 10 ' 11 : 14 and also 
the particular geometry of the het erostructure & 1 1 ' 14 

The unconventional proximity effect at SSC-FM 
interfaces clearly evidences an intimate connection 
between magnetism and spin-triplet pairing. Such 
devices can provide only a limited understanding of 
this interplay, however, as the FM determines the 
triplet pairing correlations. To be able to control 
the spin triplet pairing independently, it would be 
necessary to replace the SSC with a triplet supercon- 
ductor (TSC). Since the discovery of triplet super- 
conductivity in Sr 2 Ru04 ) 15 : 16 there has been steadily 
growing interest in the properties of TSC heterostruc- 
tures . 17 - 18 - 19 - 20 - 21 - 22 - 23 ' 24 - 25 - 26 - 27 ' 28 - 29 - 30 - 31 - 32 - 33 - 34 - 35 

Despite the likely intimate connection between the two 
phases, the study of devices combining TSCs and FMs is 
still in its infancy^o^i^iiH^i 34 ^ Even so, several 
exotic effects have already been predicted, such as a 0-7T 
transition in a TSC-FM-TSC Josephson junction caused 
by the mis-alignment of the vector order parameters 
(the so-called d-vectors) of the TSCs with the moment 
M of the FM tunneling barrier. 25 : 32 ^ 

The origin of this unconventional behaviour is the cou- 
pling of the FM moment to the spin of the triplet Cooper 
pair.' 55 More generally, the extra degree of freedom pro- 



vided by the Cooper pair spin is responsible for novel 
spin transport properties of TSC heterostructures. For 
example, a number of authors have demonstrated that 
a Josephson spin current flows between two TSCs when 
their d vectors are misaligned ! 23 ' 24 : 27 It has recently been 
established that a spin current may also be produced by 
the inclusion of a FM tunneling barrier in a TSC Joseph- 
son junctioni 32 : 33 : 34 : 35 Two basic mechanisms have been 
identified: the barrier can act as a spin-filter, preferen- 
tially allowing the tunneling of one spin species of Cooper 
pair over the other; alternatively, the barrier moment 
can flip the spin of a tunneling Cooper pair, which then 
acquires an extra spin-dependent phase. However pro- 
duced, the tunneling spin currents are always dependent 
upon the phase difference between the TSC condensates 
on either side of the junction. 

A bulk phase-independent contribution to the spin cur- 
rent in a TSC-FM-TSC Josephson junction was predicted 
in Ref. [35J, and subsequently also identified in a SSC-FM 
heterostructure in Ref. [T3 . The origin of this spin current 
was shown to be the spin-dependent phase shift acquired 
by the flipping of a triplet Cooper pair's spin upon re- 
flection at a (5-function-thin FM barrier. A number of 
properties were deduced: the spin current is polarized 
along the direction d x M; the current in the tunneling 
limit is (x sin(2a), where a is the angle between d and 
M; and the sign of the current displays a pronounced de- 
pendence upon the orbital structure of the bulk TSC, due 
to the orbital-dependent phase shift experienced by the 
reflected Cooper pairs. Most remarkable is that although 
this spin current is carried by Cooper pairs, and hence is 
indistinguishable from the tunneling Josephson spin cur- 
rent, it is independent of the material on the other side 
of the FM barrier. 

It is a natural question to ask if this effect also occurs 
at the interface between a bulk TSC and FM. In study- 
ing such interfaces, most authors have only addressed 
the case d || M when the spin current described above 
is not expected to occuri 20 ' 30 : 31 Although Hirai and co- 
workers considered arbitrary orientation of d and M in 
their study of the TSC-FM interface, they did not ex- 
amine the spin transport properties of the device^ It 
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FIG. 1: (color online) Schematic representation of the device 
studied in this work. We consider a junction between a bulk 
triplet superconductor and a bulk ferromagnet, separated by 
a thin insulating tunneling barrier. The d-vector of the TSC 
defines the x-axis, while the moment M of the FM lies in the 
x-y plane at an angle a to the x-axis. 



is therefore the purpose of this paper to examine the 
occurrence of a spontaneous spin current in a TSC-FM 
junction for arbitrary alignment of the TSC and FM vec- 
tor order parameters. Using a quasiclassical technique, 
we demonstrate the existence of a spin current with the 
same dependence upon the bulk TSC orbital structure 
and the relative misalignment of d and M as predicted 
in Ref. l3ol The spin current is also found to display 
a strong dependence upon the exchange-splitting of the 
FM, which is explained due to the angular-dependence of 
the spin-flip reflection probability. We discuss similarities 
and differences to the usual Josephson effect, in partic- 
ular arguing that low-temperature anomalies in the spin 
current imply a role for zero energy states at the interface 
in the transport. 



II. THEORETICAL FORMULATION 

A schematic diagram of the device studied here is 
shown in Fig. It consists of a bulk TSC and FM, sep- 
arated by a thin insulating barrier at z = 0, which we ap- 
proximate by a delta function of height U. Both materi- 
als are assumed to be in the clean limit. The Bogoliubov- 
de Gennes (BdG) equation describing the quasiparticle 
states with energy E is written in Nambu-spin space as 

fMr) A 41W) = £*M (1) 



VAt(r) -H£(t) 



where the caret indicates a 2 x 2 matrix in spin-space. 
The non-interacting Hamiltonian is 



Ho(r) 



h 2 v 2 

2m 



U6(z) 



In the interests of simplicity, we will assume that the ef- 
fective mass m is the same in the TSC and the FM>^ 
The moment of the FM is parameterized as M = 
M [sin(/3) cos(a)e x + sin(/3) sin(a)e y + cos(/3)e 2 ]. Since 
we do not include spin-orbit coupling, we can assume 
without loss of generality that (3 = 7r/2, i.e. the moment 
lies in the x-y plane. Other spin-orientations can be ob- 
tained by appropriate rotation of the system about the 
x-axis in spin space, which leaves the spin state of the 
TSC unchanged. 

The gap matrix in Eq. ([T]) is A(r) = i[a ■ d(r)]<T y where 
d(r) is the vector order parameter of the TSC. We will re- 
strict ourselves here to equal-spin-pairing unitary states, 
for which d(r) = A(r)8(— z)e x is a suitable choice. In 
such a state, the triplet Cooper pairs have ^-component 
of spin S z — ±h but the condensate has no net spin. 
The magnitude of the gap is assumed to be constant 
throughout the TSC. Due to the triplet spin state, the 
gap must reverse sign across the Fermi surface, imply- 
ing an odd-parity orbital wavefunction. We will mostly 
be concerned with three different orbital pairing states: 
Py-wave, A k = A(T)k y /k F ; p 2 -wave, A k = A(T)k z /k F ; 
and a chiralp z -l-ip^-wave state A k = A(T)[k z + ik y ]/k F . 
These different gaps are illustrated in Fig. ((2|). The gap 
magnitude A(T) displays weak-coupling temperature de- 
pendence, with T = value Ao- 

Both the TSC and FM are assumed to have circu- 
lar Fermi surfaces. For the sake of clarity, we assume 
the TSC and FM to have identical Fermi energies E F £& 
In the TSC, the Fermi surface in the normal state is 
spin-degenerate with radius k F — y // 2mE F /h 2 . If the 
exchange-splitting of the bands is less than the Fermi en- 
ergy in the FM, we have a majority spin (aligned parallel 
to M, s = +) and a minority spin (aligned anti-parallel 
to M, s = — ) Fermi surface of radius 



2 m 



-r?(E F + sgn B M) = kpVl + sX (3) 
nr 



where 



A = 9 hbM/Ef 



(4) 



is the ratio of the exchange-splitting to the Fermi en- 
ergy. The minority spin Fermi surface disappears when 
the exchange-splitting exceeds Ep , i.e. the FM is a half- 
metal. 

Solving the BdG equations Eq. ([T]), we construct the 
wavefunction for a spin-er electron-like quasiparticle with 
wavevector k incident upon the FM layer from the TSC 
using the Ansatz 

= e(-*)|*£*>e*-'+ £ [w*^e*-* 



-b >$ T ? C e- lkr 



]} 



i-.9^ B <T-Me(z) (2) 



+9(z) ]T {c eCTS <f e-*'" r + 4 ffs <&™e**"}(5) 
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FIG. 2: (color online) Schematic representation of the wave- 
function $ e<r for a spin-<r electron-like quasiparticle with 
wavevector k incident upon the FM from the TSC. Note the 
different lengths of the wavevectors in the FM due to the 
exchange-splitting of the Fermi surface. In the bottom left 
corner, we show the orbital arrangement of our three choices 
of pairing symmetry. In the p y and p z cases, the black and 
white lobes of the p-wave orbital indicate opposite signs; for 
the pz + ip y state, the gap magnitude is constant and the 
arrow indicates the direction of increasing phase. 



A schematic representation of the wavefunction \& ecr is 
shown in Fig. @. For z < 0, the Ansatz Eq. ([5]) describes 
an Andreev-reflected spin-cr' hole-like quasiparticle with 
wavevector k and reflection probability amplitude a eaa > , 
and a spin-cr' electron-like quasiparticle undergoing spec- 
ular reflection with wavevector — k = (k x , k y , —k z ) and 
probability amplitude b eatT i . Note that we make the 
standard "quasiclassical" assumption that E <C Ep, and 
so the magnitude of the wavevectors for the electron- 
like and hole-like quasiparticles are approximated to be 
identica h 37 ' 38 For z > 0, the transmission probabil- 
ity amplitudes for hole and electron quasiparticles with 
spin projection s — ± along the direction of the fer- 
romagnetic moment are c eas and d eas respectively. As 
shown in Fig. @, due to the spin-splitting of the FM 
Fermi surface, the trajectories of transmitted s — + and 



s = — quasiparticles are not coincident. Since transla- 
tional invariancc is satisfied along the x and y directions, 
the component of the wavevector parallel to the inter- 
face is preserved during the scattering. We hence have 
kps'm(0) = kps'm(6 s ). For sin(6') < k F /kp propagating 
solutions for both spin polarizations exist, and the z- 
component of k s and k s is then k SjZ — kpy // cos 2 (8) + sX. 
Since kp < kp < kp, there is however a critical an- 
gle 9 C = arcsin(fc^/fc_F) such that for 6 C < \9\ < n/2 
the z-component of the wavevector of the transmitted 
s = — quasiparticles is purely imaginary. The result- 
ing evanescent wave is exponentially suppressed on an 
inverse length scale k — kp^X — cos 2 (6). 

The wavefunction Ansatz Eq. ([5|) is expressed in terms 
of the spinors for the bulk TSC and FM phases. In the 
FM we have the spinors for electrons and holes with spin 
s = ± 



FM 



(se-"7%/2,l/V2,0,0 



<f = (O, 0,^/^,1/^2 



(G) 
(7) 



For the TSC, the spinors $ke(/s) a ^ or an electron-like 
(hole-like) quasiparticle with spin a and wavevector k 



foTSC 

<f> TSC 
<f> TSC 

qTSC 



(s k u k ,0,-u k ,0) T 
(skWk,0, -M k , 0) T 
(0,s k u k ,0,-u k ) T 
(0, s k v k ,0, u k ) T 



(8) 
(9) 
(10) 

(11) 



where u k = y/(E + Sl k )/2E, w k = y/(E-Q k )/2E, ft k = 
I 2 and s k = A k /|A k |. 

are determined 



V E2 - l A ki 

The probability amplitudes in Eq. 
by the boundary conditions obeyed by the wavefunction 
at the TSC-FM interface. In particular, we require that 
the wavefunction is continuous at the interface, i.e. 



^e<r(r)| z= o- = * e <r(r)| 



z=0+ 



(12) 



and also that the first derivative of the wavefunction 
obeys the condition 



c"F eCT (r) 



dz 



=0+ 



dz 



Zfc F * eCT (r)| z=0 (13) 



where Z — 2mU/h 2 kp is a dimensionless constant char- 
acterizing the height of the insulating barrier. These con- 
ditions yield eight coupled equations for the probability 
amplitudes. In general, the probability amplitudes are 
functions of the incident wavevector k and the quasipar- 
ticle energy E. 

It is straight-forward to modify the Ansatz Eq. ([5]) for 
the wavefunctions ^^(r) describing the scattering of a 
hole-like quasiparticle incident upon the FM from the 



TSC. The probability amplitudes aha 



etc. for 



4 



this case are determined by applying the same boundary 
conditions as for ^ e(T (r). 

The currents in the TSC may be evaluated using the 
generalization of the Furusaki-Tsukuda formula to triplet 
pairing i 24 i 39 For the currents flowing perpendicular to the 
interface, we hence express the charge current Ic and the 
z- and y-components of the spin current, Is, z and Is, v 
respectively, in terms of the Andreev reflection probabil- 
ity amplitudes a e th)crcr' with energy argument analytically 
continued to iui n ' 



T=0AT 



r=0.004r 



Ic 



e 

2Ti 



( dk z dk v Q(k z )^ VV 




0Wr(k, iu n ) 



n 



n, — k 



IS,: 



Is,y 




-a haa Ck,iu! n ) >(15) 



-k 



dk z dk y (d(k z )j^-— a 

n a 



a e<jS (k, iu> n ) 



k F (3h 
IA 



-a hcr5 .(k,iw„) >(16) 



where ui n = (2n - l)n /f3, f2„,k = v /u} n + l^k| 2 and 
a = —a. The expression for the x-component of the 
spin current is vanishing, which reflects the fact that the 
Cooper pairs in the TSC do not have a spin component 
parallel to the d-vector^ As required for the TSC-FM 
junction, the charge current is found to vanish in all con- 
sidered circumstances. Furthermore, the y-component of 
the spin current is also vanishing for the choice of M 
adopted here. As mentioned above, any orientation of 
M can be rotated into the x-y plane without changing 
the pairing state of the TSC. Since the polarization of the 
spin current must also be correspondingly rotated under 
such a transformation, we hence conclude that the spin 
current is always polarized along the direction d x M. 35 

Due to the time-reversal-symmetry-breaking pairing 
state in the p z + ip y -~w&ve TSC, here we expect to find 
surface charge and spin currents flowing parallel to the 
interface. 17 ! 20 ! 34 We will not be concerned with such cur- 
rents in what follows, presenting only results for the spin 
current flowing perpendicular to the interface. For con- 
venience, we shall also adopt units where h = 1. In all 
plots we take the interface parameter Z = 1. 
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FIG. 3: (color online) Spin current as a function of the angle 
a for the p y (black solid curve), p z (red dashed curve), and 
Pz + ipy-wave (blue dot-dash curve) TSCs for various values 
of the exchange-splitting A. The left column shows the spin 
currents at T = 0.4T C while the right column shows the spin 
currents for T = 0.004T C . 



III. RESULTS 

We find that the spin current is an odd periodic func- 
tion of a with period n. As shown in Fig. ((3|), for most 
situations the spin current obeys Is, z oc sin(2a). This 
dependence upon a is reminiscent of the usual Joseph- 
son charge current vs phase relationship for Josephson 
junctions, suggesting that the angle of misalignment be- 
tween d and M plays a role similar to the phase differ- 
ence. As was argued in Ref. HH, this analogy is valid: 
a spin a Cooper pair incident upon the barrier acquires 
a phase shift — 2aa when undergoing a spin-flip at re- 
flection. Interpreting spin-flip reflection as "tunneling" 
between the spin-t and spin- J, condensates of the TSC, 
the phase shift — 2a a is therefore the effective phase dif- 
ference for a Cooper pair "tunneling" between the spin a 
and spin —a condensates. Naturally, this drives a Joseph- 
son charge current of equal magnitude but opposite sign 
in each spin sector of the TSC, thus producing a finite 
spin current but vanishing total charge current. Note, 
however, that the phase difference in a Josephson junc- 
tion is a property of the wavefunctions of the bulk con- 
densates on either side of the tunneling barrier; in the 
TSC-FM junction, in contrast, the effective "phase dif- 
ference" results entirely from a property of the interface. 

Higher harmonics in 2a clearly appear in the spin cur- 
rent for the p z -wave case at low temperatures and moder- 
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FIG. 4: (color online) (a) Definition of the angle 7 parame- 
terizing the alignment of the p-wave orbital in the TSC with 
the interface. 7 = gives the p^-wave orbital state, while 
7 = 7r/2 gives the p z -wave orbital state, (b) Spin current as 
a function of 7 for various values of A. We take a = O.ln and 
T = 0.4T C . (c) The spin current as a function of the incident 
angle 6 in the p z + ip y junction for a = 0.17T and T = 0.4T C ; 
(d) shows the same curves at T = 0.004T C . 
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FIG. 5: (color online) (a) Spin current as a function of the 
ratio A = q^bM/Ef. Note that the current for the p z + ip y - 
wave TSC changes sign at A ~ 0.4. (b) Spin current at low 
A. Curves have identical meaning as in panel (a). In both 
(a) and (b) we set a = O.lvr and T = 0AT c . (c) T > T c 
normal state spin-slip reflection probability as a function of 
the incident quasiparticle angle 6. 



ate to strong exchange splitting A > 0.5. Continuing the 
Josephson junction analogy, we interpret this as evidence 
of coherent spin-flip reflection of multiple Cooper pairs. 
In a Josephson junction, such processes are associated 
with the formation of a zero energy bound state at the 
tunneling barrier, which allows the resonant tunneling of 
multiple Cooper pairs. It is well known that zero energy 
states form at unconventional superconductor interfaces 
when the orbitals are aligned such that a quasiparticle 
specularly reflected at the interface experiences a sign 
reversal of the superconducting order parameter^ Such 
a state does not therefore occur in the p y junction, while 
it is present for all k in the p z junction, and present only 
for k = ke z in the p z + ip y junction. This is consistent 
with the absence of higher harmonics in the spin current 
for the p y and p z + ip y junctions. Although suggesting 
an important role for the zero energy states in the spin 
current, it does not explain the absence of the higher 
harmonics at low A. This indicates that the coherent re- 
flection of multiple Cooper pairs is also controlled by the 
bulk properties of the FM. 

It is clear from Fig. ([3]) that the spin current is strongly 
influenced by the orbital structure of the TSC; in par- 
ticular, the spin currents in the p z and p y junctions al- 
ways have opposite sign. In the former case, the reflected 
Cooper pairs acquire an additional it phase shift due to 
the reversed sign of the superconducting gap, thus re- 
versing the spin current relative to the p y case. This 
can be strikingly demonstrated by examining the varia- 
tion of the spin current as the p-wave orbital in a time- 
reversal-symmetric TSC is rotated from the p y -wave to 
the p z -wave configurations. We parameterize the orien- 



tation of the orbital by the angle 7 between the p-wave 
orbital maximum and the y-axis, see Fig. ©(a). Only 
Cooper pairs with incident angle < 6 < 7 experience a 
reversal of the gap sign upon reflection, and hence con- 
tribute a spin current of opposite sign relative to the 
Cooper pairs with incident angle 7 < 6 < ir/2. As shown 
in Fig. ((4]) (b) , the former contribution to the spin current 
comes to dominate the latter as 7 is increased from to 
7r/2, with the spin current consequently changing sign at 
some critical value of 7. 

A similar effect occurs in the p z + ip y case, but here 
the additional orbital phase shift experienced by the re- 
flected Cooper pairs is tt — 2 arctan(/cj,/fc z ), i.e. it de- 
pends upon the incident trajectory. Normally incident 
Cooper pairs thus undergo an additional tt phase shift 
as in the p z junction, whereas a trajectory grazing the 
interface has no additional phase shift as in the p y junc- 
tion. The angle of incidence therefore determines the sign 
of the spin current contributed by the reflected Cooper 
pair. Making the change of variables k z — > kpcoa(9), 
k y — > fci?sin(#) in Eq. (|15p. we define the angle-resolved 
spin current Is,z(fi) by 

,tt/2 

Is,z= / dOIsAO) (17) 

J-tt/2 

We plot Is,z(8) in Fig. dl])(c) and (d), where we see that 
it changes sign at 9 » 0.2ir. Note that for A < I the 
angle-resolved spin current is sharply peaked at the criti- 
cal incident angle C , where the transmitted s = — quasi- 
particle has vanishing z-component of its wavevector. 
The variation of the spin current with A is shown 



in Fig. ©(a). There are some common features for all 
choices of the TSC orbital: at low A the spin current 
goes as ~ A 2 ; a maximum value of the spin current is 
reached at (p z and p z + ip y ) or close to (p y ) A = 1, 
where the minority spin Fermi surface disappears; and 
for A > 1 the magnitude of the spin current displays slow 
monotonic decrease. Despite these similarities, the three 
curves nevertheless have some key distinguishing features 
at A < 1 . In the p y case, the spin current increases much 
more slowly with A than for the other junctions, while 
the maximum is smooth and occurs before the disap- 
pearance of the minority Fermi surface. Although the 
spin currents in the p z and p z + ip y junctions show sim- 
ilar A-dependence for A > 0.5, the spin current in the 
Pz+ipy junction changes sign at A w 0.39 [see Fig. ©(b)]. 
We further note that the spin current in the p y junction 
has greater magnitude than that in the p z junction for 
A < 0.2, but the smallest magnitude of all the currents 
for A > 0.65. 

The transport properties of a junction can be formu- 
lated in terms of the normal state scattering matrix?! 
Insight into the complicated dependence of the spin cur- 
rent upon A can therefore be gained by examining the 
T > T c spin- flip reflection probability |i? s /| 2 , given by 
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FIG. 6: (color online) The maximum value of the spin current 
obtained by varying a for (a) A = 0.05, (b) A = 0.5, (c) 
A = 1.0 and (d) A = 1.5. The curves in panels (b)-(d) have 
the same meaning as those in panel (a); note that the p z +ip y 
curve is multiplied by 5 in (b) for clarity. 
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fc^(fc+,z-fc-,z) 2 101 < n 

fc;K,+* a ) 

[fc2 + (Zfc F +K)2][Z2fcJ. + (fcz+fc+,z) 2 J 

(18) 

Note that |i? s /| 2 is finite when Z = 0, and so the spin cur- 
rent survives when the insulating barrier is removed. We 
plot |i? s/ | 2 for Z = 1 in Fig. ©(c). When A < 1 the spin- 
flip reflection probability is sharply peaked at the critical 
angle 6 C , and grows in magnitude as A is increased. For 
A > 1 the reflection probability remains peaked at 6 = 0, 
although it decreases from the A = 1 maximum. For 
larger values of Z, |-R S /| 2 may very slightly increase with 
increasing A at 9 ^ (not shown). Comparison of the 
angle-resolved spin current Is,z(9) in Fig. ©(c) and (d) 
with Fig. ©(c) clearly shows the connection of |i? s /| 2 
to the transport: the spin-flip reflection of Cooper pairs 
with incident angle 9 close to 9 C is strongly favored, and 
tends to dominate the spin current. The sign change of 
the spin current in the p z + ip y junction with increasing 
A can thus be understood as arising from the shift in the 
peak of \R s f\ 2 ' at A <C 1 this strongly favors grazing 
trajectories experiencing only small orbital phase shift, 
while at A ~ 1 almost-normal trajectories with orbital 
phase shift close to n dominate. 

The gap anisotropy in the p z and p y junctions implies 
a greater contribution to the spin current for trajectories 
along which the gap magnitude is maximal. The peak in 
|i? s /| 2 at grazing trajectories for A <C 1 therefore results 
in a greater spin current in the p y junction than the p z 
junction. As A is increased, the reduction of 9 C leads 
to the observed strong enhancement of the spin current 
in the p z junction. Because the reflection probability at 
high 9 also grows with A, the spin current continues to 



increase in the p y junction, albeit at a much slower rate 
than in the p z junction. The reduction of the spin current 
at A > 1 is due to the reduction of the peak in |i? s /| 2 at 
6» = Q. 

The last quantity of interest is the maximum spin cur- 
rent max{/s iZ } obtained by varying a. This is plotted as 
a function of T for fixed A in Fig. © . As can be seen, for 
all A the magnitude of the spin current in the p y junction 
shows a moderate enhancement with decreasing temper- 
ature, with an apparent plateauing at T < 0.1T C . In the 
p z junction, in contrast, there is a strong enhancement 
of the spin current magnitude with decreasing tempera- 
ture. The temperature dependence of max{is iZ } in the 
p z and p y cases is reminiscent of the critical current in 
Josephson junctions with and without zero energy states 
respectively ) 21 ! 38 again indicating a role for such states in 
the device studied here. The temperature dependence of 
max{/,g,z} in the p z + ip y junction is determined by A: 
at low and high A wc find the moderate increase with de- 
creasing T characteristic of the p y junction, and also seen 
in the critical current of non-magnetic Josephson junc- 
tions between p z + ip y TSCsJ^ At intermediate A, how- 
ever, we find a non-monotonic dependence of max{7s,z} 
upon T . This can be understood by examining the angle- 
resolved spin current in Fig. ©(c) and (d): for A = 0.5, 
we observe a particularly strong enhancement of the cur- 
rent for trajectories near 9 C as the temperature is lowered. 
Since the current for such trajectories has opposite sign 
to that of the total current, there is consequently a re- 
duction of the total spin current. For a small range of 
A ~ 0.4, the spin current in the p z +ip y junction reverses 
sign with decreasing temperature. 
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FIG. 7: (color online) The spin-flip reflection of the Cooper 
pairs imparts the spin change SS || d x M to the FM. We 
speculate that this results in a precession of the magnetic 
moment around the axis defined by the d vector. 

IV. CONCLUSIONS 

This work presents an analysis of the spontaneous spin 
current generated in a TSC by contact with a bulk metal- 
lic FM. The essential requirement for the appearance 
of this spin current is that the vector order parameters 
of the TSC and FM are misaligned but not mutually 
perpendicular. Following Ref. l35l . the spin current has 
been interpreted as a Josephson-like effect, where the 
phase shift picked up by Cooper pairs undergoing spin- 
flip reflection at the interface plays the role of the phase 
difference between the superconductors in a Josephson 
junction. This analogy is supported by the temperature 
evolution of the maximum spin current, and the low- 
temperature dependence upon a, which suggests that 
resonant reflection of multiple Cooper pairs through zero 
energy states plays a significant role in the current. 

The spin current nevertheless possesses several proper- 
ties which are not anticipated by the Josephson junction 
analogy. As it arises from reflection processes, the choice 
of orbital pairing state in the TSC can determine the 
sign of the spin current, due to the orbital phase shift 
experienced by the reflected Cooper pairs. In the case of 
the p z + ip y junction, this produces the interesting result 
that spin-flip reflected Cooper pairs with different inci- 
dent trajectories can carry spin currents of opposite sign. 
The exchange-splitting of the FM also influences the spin 
current, most significantly when the FM has both a mi- 
nority and majority spin Fermi surface. This is closely 
connected to the orbital structure of the TSC, as the 
exchange-splitting determines the angular-dependence of 
the spin-flip reflection probability. 

A few critical remarks upon our method are necessary. 
Considering that interface effects play the essential role 
in the generation of the spin current, it is reasonable to 
question our approximation that the TSC and FM order 
parameters are both constant up to the insulating bar- 
rier. It is well known that there is a strong suppression of 
the p z -w&ve gap at partially-transparent interfaces ; 17 ! 26 



but as this does not change the spin structure of the 
TSC it is only likely to renormalize the results presented 
here. More interesting is a possible change in the orien- 
tation of d(r) and M(r) close to the interface due to an 
unconventional proximity effect. Although the analysis 
of Ref. HU indicates that the spin current will survive a 
self-consistent treatment of the superconductor, it is de- 
sirable that the TSC and FM be treated on an equal level. 
This is possible, for example, within a real-space Hartree- 
Fock analysis^! 2 ^. The TSC-FM interface has already 
been studied using this method ) 20 ' 31 but only for the case 
d || M which we predict to display vanishing spin current. 
Much scope therefore remains for further investigations. 

Another important assumption is that the magnetic 
moment of the FM is constant. As pointed out in Ref. l35l. 
however, the spin-flip reflection of the Cooper pair im- 
parts a spin <5S || d x M to the FM. If \SS\ < |M|, the 
effect of the imparted spin should be to slightly rotate 
the magnetic moment about the axis defined by d, with- 
out significantly altering its magnitude, see the cartoon 
representation in Fig. (0. This is a generalization of the 
well-known spin-transfer torque effect of spintronics to a 
superconducting device4°- For a FM without anisotropy, 
we speculate that the cumulative effect of many spin-flip 
reflection events is to cause the precession of the mag- 
netic moment. If this precession is sufficiently slow, the 
results presented here should still remain valid, and we 
anticipate a periodic modulation of the polarization of 
the spin current due to the rotating moment. Rigorous 
verification of these claims requires a more sophisticated 
analysis than our argument above, and is left for future 
work. 

Finally, we consider the prospects for the experimen- 
tal verification of our predictions. Direct measurement of 
spin currents in superconductors is unlikely to be easy, al- 
though several proposals exist, e.g. spin-resolved neutron 
scattering 41 , or ARPES with circularly-polarized light. 42 
Alternatively, one could search for evidence of a spin- 
accumulation at the edge of the TSC opposite to the 
interface with the FM. If our speculation above proves to 
be well-founded, it might also prove possible to deduce 
the existence of a spin current by measuring the preces- 
sion of M. Initial characterization of a TSC-FM device 
would, however, almost certainly involve measurements 
of the LDOS at the interface 3 ^ and the tunneling con- 
ductance^ The unique behaviour of these quantities at 
TSC-FM interfaces is a signature of the unconventional 
interplay of ferromagnetism and triplet superconductiv- 
ity, of which the spontaneous spin current is only one 
manifestation. 
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